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We numerically study two-dimensional quantum turbulence with a Gross-Pitaevskii model. With 
the energy initially accumulated at large scale, quantum turbulence with many quantized vortex 
points is generated. Due to the lack of enstrophy conservation in this model, direct energy cascade 
with a Kolmogorov-Obukhov energy spectrum E{k) oc k~^^^ is observed, which is quite different 
from two-dimensional incompressible classical turbulence in the decaying case. A positive value for 
the energy flux guarantees a direct energy cascade in the inertial range (from large to small scales). 
After almost all the energy at the large scale cascades to the small scale, the compressible kinetic 
energy realizes the thermodynamic equilibrium state without quantized vortices. 

PACS numbers: 67.25.dk,47.37.-|-q 
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INTRODUCTION 



The experimental discovery of Bose-Einstein conden- 
sates (BEC) 15 years ago [1-3], long after their theo- 
retical prediction in 1924-25 [4, 5], has renewed interest 
in this field. BEC systems are of great interest as they 
promise the opportunity to study new nonlinear dynam- 
ical systems built with a high degree of control and flex- 
ibility. Also, theoretical and numerical studies of BEC 
systems are of general importance for nonlinear physics, 
as BEC systems are described by one of the most im- 
portant and universal partial differential equations, the 
nonlinear Schrodinger equation, called in this field the 
Gross-Pitaevskii (GP) equation (GPE) [6, 7]: 
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v-^* - -.g|*r* = 



(1) 



Here the condensate wave function ^(r,t) plays the role 
of the complex order parameter. 

The GPE (1) describes a Bose gas at low tempera- 
ture, which may behave similarly to a superfluid, and 
can describe its random (turbulent) motion, i.e. quan- 
tum turbulence (QT). QT physics, comprising tangled 
quantized vortices, is an important research topic in low- 
temperature physics [8, 9]. Stimulated by recent experi- 
ments on both superfluid '^He and superfluid ^He, where 
a few similarities have been observed between quantum 
and classical turbulence [10-17], studies on QT have en- 
tered a new stage where one of the main motivations is to 
investigate the relationship between quantum and classi- 
cal turbulence. In particular, the Kolmogorov-Obukhov 
turbulent kinetic energy spectrum -E(A:) cx k~^l'^ has 
been observed in laboratory experiments on superfluid 
'^He similar to that in normal fluids [18]. The phys- 
ical explanation for this is very simple: for superfluid 
motion at scales £ essentially exceeding the mean inter- 
vortex distance I, the quantization of the vortex lines 
can be neglected. This range of scales corresponds to 
the quasi-classical limit and this type of turbulence can 



be called quasi- classical turbulence (QCT). QCT in the 
three-dimensional (3D) case has subsequently been found 
in numerical simulations with the GPE as well as in the 
vortex filament model [19] that can describe motions of 
quantized vortex lines with a prescribed core structure. 

It is well known that Euler equations, that describe 
motion of ideal (inviscid) fluids, in the two-dimensional 
(2D) case exhibit an additional second quadratic motion 
invariant, enstrophy fl, as well as kinetic energy E: 



n= u'^dr, uj = [V X v]^ . (2) 



E^- I \v\^dr, 



Kraichnan recognized [20] that the appearance of 
drastically modifies the physics of 2D turbulence. It 
changes the direction of the energy flux in the classi- 
cal Richardson-Kolmogorov cascade. Instead of a direct 
energy cascade (from large to small scales), like in 3D 
turbulence, in the 2D case there is an inverse energy 
cascade, from small to large scales. In both the 3D- and 
2D-cases the energy spectrum is the same, E(k) oc k~^/^. 
In 2D turbulence the direct energy cascade is replaced by 
the direct enstrophy cascade with the energy spectrum 
E{k) cx k^^ (subject to certain logarithmic corrections, 
unimportant in the current discussion). These predic- 
tions have been confirmed in laboratory experiments, see 
e.g. Ref. [21], and large scale direct numerical simulations 
(DNS) of the Navier-Stokes equations, see e.g. Refs [22- 
24]. 

An important issue is whether the inverse energy cas- 
cade is a feature of 2D QT. The qualitative answer given 
in our paper is "not necessarily" . The physical reason for 
this is quite simple and general. 

In superfluids, the enstrophy coincides (up to a pref- 
actor) with the total number of quantized vortex points 
[25] . In GP dynamics the total number of vortices is not 
conserved. They can appear in pair-creation or disap- 
pear in pair-annihilation. Therefore, Kraichnan's argu- 
ments [20] become, generally speaking irrelevant for 2D 
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QT. Thus, for a certain range of parameters, where the 
creation/annihilation of vortex pairs become dynamicahy 
important, we can expect a direct energy cascade in 2D 
QT, exactly as in 3D classical turbulence. This is an 
important conclusion of our paper. Simulating the free 
evolution of 2D GPE (from an initial condition with en- 
ergy located at large scales) we observed a direct energy 
cascade with the Kolmogorov-Obukhov energy spectrum 
E{k) oc 

The same argument can be made based on Euler equa- 
tions. Indeed, the superfluid density in the vicinity of 
a vortex core is different from that in a vortex-free su- 
perfluid. Therefore the creation/annihilation of vortex 
pairs lead to variations of the superfluid density. In other 
words, the GPE is reduced to a compressible Euler equa- 
tion. However, a compressible Euler equation does not 
preserve enstrophy. Therefore at some level of compress- 
ibility (characterized by the Mach number, the ratio of 
the turbulent velocity fluctuations to the sound velocity) 
the direction of the energy flux can change sign and in- 
stead of an inverse energy cascade, we observe a direct 
cascade, typical for 3D turbulence. 

In this study, we observe the direct energy cascade 
in numerical experiments with 2D GPE. Starting from 
initial conditions with energy localized at large scales 
we observed intermediate asymptotic Richardson cas- 
cade transporting energy of the system toward small 
scales. The energy spectrum at this stage is close to 
the Kolmogorov-Obukhov law E{k) oc k~^/^. At later 
stages, when an essential part of the energy reaches the 
smallest scales available in our simulations (mean inter- 
vortex spacing, which is a few times larger than the core 
diameter), the system quickly evolves toward thermody- 
namic equilibrium. As an independent test of the direc- 
tion of the energy flux we analyzed the energy balance 
equation in fc-spacc and found that the sign of the energy 
flux indeed corresponds to tlw^ dircK;t cascade. 

In addition, we studied the properties of the finite 
steady-state of the system by a different means, e.g. by 
analyzing the power spectra in the frequency domain for 
motion with different wave vectors. We show that the 
position of the maximum of the power spectra a;max(fc) 
depends on the dispersion relation derived from Bogoli- 
ubov's microscopic theory [6, 7]. 

In addition, we investigate how the time evolution of 
the system depends on the level of nonlinearity by quali- 
tatively comparing its behavior with the same initial con- 
ditions but different coupling constants g. 

The paper is organized as follows. After this introduc- 
tion, in Sec. II we explain the analytical and numerical 
formulations of the GP model of compressible quantum 
turbulence. Next, in Sec. Ill we explain the turbulent 
state, i.e. the time development of the energy, kinetic 
energy spectra, and incompressible kinetic energy flux. 
We also show how the direct energy cascade with the 
Kolmogorov-Obukhov law is formed and prove the direct 
energy cascade with the use of the approximated incom- 
pressible kinetic energy flux. Then, through numerical 



analysis of the power spectrum of the compressible effec- 
tive velocity, we find that Bogoliubov's microscopic the- 
ory holds in the thermodynamic equilibrium state full of 
compressible kinetic energy. Finally, in Sec IV we state 
the conclusions of the paper. 



II. GP MODEL OF COMPRESSIBLE 
TURBULENCE 

A. GP and "quantum" Euler equations 

It is well known [6, 7] that GPE conserves the total 
energy (Hamiltonian) H and "particle number" N: 



H = 
N = 



|V*| 



-MM 



dr , 



dr . 



(3a) 
(3b) 



By analogy with the Schrodinger equation we can con- 
sider the probability |^'(r,t)|^ as a particle density: 



pir,t)^\<SJir,t)\'' . 



(4a) 



Then the conservation law of N can be considered as 
the conservation of the total mass of the system M = 
J pdr = N (assuming that the "particle mass" is unity), 
which can be written as a continuity equation: 



dp „ . 



. 



Here j is the particle flux 



3 = -(*V** - **V*) , 



(4b) 



(4c) 



which can be presented in the familiar form j = pv with 
the "fluid" velocity: 



v = ve . 



(4d) 



Here the phase 9 is deflned via the Madelung transfor- 
mation: 



* = y/pexp{ie), 



(4e) 



which maps the GPE (1) to the Euler equation for an 
ideal compressible fluid with an extra quantum pressure 
term. 



B. Decomposition of the system energy 

Following Ref. [16] it is instructive to decompose the 
total energy of the system H, Eq. (3a), conserved by 
GPE (1), into four parts. The first part originates 
from the interaction energy in the Hamiltonian (3a) of 
the GPE (1) with the density 5(|*|^/4 = gp'^/A. Bear- 
ing in mind that J pdr =const., we can introduce the 
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density of the "internal energy" as gp^/A counted from 
(5/2) (p-i): 



(5) 



Next, we define the density of the fluid kinetic energy 
as usual: 

^kin(r) = ^p|t;p = w = y/pv. (6a) 

This kinetic energy can be divided into compressible and 
incompressible parts by decomposing the effective "ve- 
locity" field w into divergent free and potential w'^ 
parts: 

^kin(r)=4inW+^&„(r), (6b) 
£Ur)^^\wf, fUr-)^^^!', (6c) 

because ■ = 0. 

Now, the first term in the Hamiltonian (3a) in the rep- 
resentation of (4d), (4e) can be presented as follows: 



2' ' 2 
^qnt(r) ^ i IVyTpf . 



(7) 



The second term here, f qnt('"), has been termed the quan- 
tum energy density [17]. 

Finally, the total (conserved) energy of the system can 
be decomposed into four parts: 

H= f [£i„t(r) + 4„(r) + £^,,{r) + S^^^ir)] dr . (8) 



C. Numerical procedure 

We solve the 2D GPE (1) by the pseudo-spectral 
method in the domain 256^ with the fourth-order Runge- 
Kutta method for the time development. For details, see 
Ref. [11]. 

To generate large-scale turbulent fiow, the initial con- 
dition of the wave function is set to the random phase 
state *(r,0) = cxp[i0(r,O)] with 



e{k, 0) 



eoex.p[ia{k)] {Ak < \k\ < 3Afc) , , . 
(otherwise) . ^ ^ 



Here 6{k, 0) = 9*{—k, 0), a{k) is randomly taken for each 
fe in the range (— tt, tt), and Afc = 2'k/L is the wavcmim- 
ber grid. The initial density is uniform, |>I'(r,0)p = 1. 
The phase distribution (9) is energetically high and un- 
stable so that many quantized vortex pairs arc created. 
With these initial conditions we obtain, during time evo- 
lution, a 2D QT composed of a random configuration of 
quantized vortices. 

To investigate the character of the 2D QT we evolve the 
GPE (1) as described above and compute the following: 



1. Time evolution of all four total energy compo- 
nents (8) per total particle number N and total 
vortex number AAqv for different values of the cou- 
pling constant g = 1 , 2 , 4, shown in Fig. 1. 

2. Time evolution of the dimcnsionlcss "level of non- 
linearity", defined as i^int / (-E'kin + -Eqnt), shown in 
Fig. 2. 

3. Compressible and incompressible kinetic energy 
spectra, -E'kin(^) -^kin(^)' shown for g = 1 , 2 , 4 
in Figs. 3, 4, and defined below: 

Kinit) = ^ J mk)\'dk^ J E-^,^ik)dk, (10a) 
EL{t) = ^ J \w\k)\'dk^ J El,^{k)dk, (10b) 
1 = -! j\^{k)\^dk = j N{k)dk, (10c) 



where 



(lOd) 



13=1,2 

L(fe)= E {s^,-^y,ik). (10 



4. The spectra E^.^{k), E'l-^ik) and N{k) for different 
moments of time and different values of g, shown 
in Figs. 3, 4, and 5 respectively. 

5. Mean incompressible kinetic energy flux in fc-space, 
s^{k) shown in Fig. 6 {g = 1 , 2, 4). 

6. Distribution of vortex-pair separations, shown for 
5 = 1 , 2 , 4 in Fig. 7. 

7. Frequency spectra of compressible and incompress- 
ible velocity components for different k and dif- 
ferent g at later time moments, shown in Figs. 8 
and 10 respectively. Positions of the maxima of the 
compressible frequency spectra Winax(^) are shown 
in Fig. 9 for different g and different moments of 
time. 



D. Incompressible kinetic energy flux 

One of our aims is to estimate from numerics the value 
and direction of the energy flux over a variety of scales in 
the incompressible and compressible subsystems, e^{k,t) 
and s'^{k,t). To do this we use the observation (see be- 
low Fig. 1) that during most of the evolution (say for 
t > 2.5) the total kinetic energy is practically conserved. 
We can then neglect the energy exchange between these 
subsystems. In contrast, there is a permanent energy 
flux from the incompressible to the compressible subsys- 
tem and this exchange has to be accounted for. Therefore 
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FIG. 1. (Color online) Time evolution of the dimensionless total energy Et, total kinetic energy -Ekin and its compressible 
iS^in S'lid incompressible parts. All energies are normalized with the particle number A'^, as shown in Eqs. (10). The time 
dependence of the total vortex number (a measure of the enstrophy) is also shown. Left panels: g = I, Middle panels: g = 2, 
Right panels: g — 4. Upper panels: initial evolution, t < 5, Middle panels: intermediate stage, t < 40, Lower panels: latest 
stage, t < 4500 for 5 = 1, i < 2000 for g = 2, and t < 600 for g = 4. 



the (approximate) balance equations for the energy den- 
sities E^^^^{k,t) and El^^{k,t), introduced by Eqs. (10), 
can be written as follows: 



dt dk 
dEl,^{k,t) , de\k,t) 



dt 



dk 



+ F"ik,t)^0, (11a) 
-F"{k,t)^0. (lib) 



Here F"{k,t) describes the energy exchange between 
subsystems. Integrating these equations we conclude: 



e%k,t) 
e'{k,t) = - 



fe 



dt 



dEUk',t) 



dt 



+ F"{k',t) dfc',(12a) 



F"{k',t) dk' .(12b) 
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A zero-order approximation is applied in the analysis of 
Eqs. (12) to eliminate F^^{k,t) in these equations. This 
approximation can be improved as follows. We know that 
energy flux terms do not contribute to the total balance: 



III. DECAY OF 2D COMPRESSIBLE 
TURBULENCE 

Free evolution of system energies and vortex 
number 



I 



dt 



dt 



+ F''\k',t) dk' = 0, (13a) 



F"(fc',t) dk' = 0, (13b) 



which allows us to find from numerics F"{k', t)dk' 

at any instant of time, for example, in the following way: 



"'max ■>'ma 

j F"{k',t)dk'= J 



^d[EL{k',t)-E^,^{k',t)] 



2dt 



femin 



femin 



(14a) 

Next we can approximate F"{k',t) in the factorized 
form 



F^\k\t)^ fit)^ik') 



(14b) 



with a prescribed fc'-dependent function (p{k') that can 
be determined by reasonable physical arguments. For 
example, assuming that F'^^{k',t) mostly originates from 
the collapse of vortex pairs, noting that the velocity 
around each vortex decays like 1/r', giving in the k'- 
representation proportional to k' and that El.-^^{k') cx 
k' \ w{k')\'^. In our calculations we took <^{k') oc k'^, or 
in the normalized form 



max 



ip{k') = 4k'yk 

Then from Eq. (14a) we find 

k 



(14c) 



Equations (14b), (14c), and (14d) can be substituted 
back into Eqs. (12) to get an improved approximation 
for the flux. Finally, we obtain the incompressible ki- 
netic energy flux in a form with (14d): 



e\k,t) 



I 



dEl,^{k',t)^^, .„,../ k 



dt 



-dk' + fit) 



(15) 



Numerical results for the energy flux, obtained in this 
way, are shown in Fig. 6 and will be discussed later. 



Energy evolutions of the GP system, starting from the 
same initial condition ^(r, 0), but with different values of 
the coupling constant, g = 1, g = 2, and (7 = 4, are shown 
in Fig. 1. Visual analysis of the figure leads to several 
qualitative statements concerning the system behavior: 

• As expected the total energy, Hamiltonian is time 
independent and has larger values for larger g (due to the 
g\^\'^ contribution). 

• At the initial stages of the time evolution (for t < 
5, upper panels) an intensive process of vortex creation 
leads to fast transformation of the compressible kinetic 
energy into incompressible kinetic energy and quantum 
energy (not shown) . The incompressible kinetic energy is 
larger than the compressible kinetic energy for the time 
interval 2<t<4. 

• At later stages {t > 5) the total kinetic energy is 
practically independent of time. Moreover its values are 
more or less the same for all g. This can be easily under- 
stood: the kinetic energy is quadratic in ^ and does not 
contain the nonlinear term proportional to g. 

• The largest maximal value of the vortex number 
Afqv{t) is achieved at a smaller value g = 1 (Left panel) 
at the crossover time t « 5. This is because for t > 5 the 
probability of the dominant nonlinear process of vortex- 
pair annihilation is larger for larger g. We see that the 
fastest decay of A/'qv(t) corresponds to 5 = 4, Right panel. 

• A faster decay of the vortex number for larger g 
leads to a larger increase in the flow compressibility with 
an increase in g. Correspondingly, the most compressible 
flow (largest value of the ratio E^-^/El^^) is reached (in 
the presented data) for g = 4, and the least for g ~ 1. 

• The system finally reaches its equilibrium state with- 
out quantized vortices for all g. It takes longer to reach 
thermodynamic equilibrium state for smaller g. 

• To clarify in more detail why a large value of the cou- 
pling g corresponds to a faster time evolution, we com- 
pute and plot in Fig. 2 the time evolution of the interac- 
tion energy i^int , defined by Eq. (5) (and proportional to 
g). This energy has to be compared with the smn of the 
kinetic and quantum energies, i^kin + -Eqnt, which origi- 
nates from the quadratic term in the Hamiltonian (3a) 
or the linear term in the GPE (1). This sum is shown in 
Fig. 2. We see that typically Eint is essentially smaller 
than E'kin + E^^t and decays in time. Therefore the ratio 



R = Eint I (-Ekin + E, 



qnt; 



(16) 



shown in Fig. 2, characterizes how far our system deviates 
from its linear approximation. Accordingly, the ratio R 
can be called the nonlinearity level. The larger i?, the 
more quickly the system should evolve to its equilibrium 
state. Indeed, we see that the maximal value of i? ~ 0.8 
is reached (at t w 0.8) for g = A, while for 5 = 1 the 
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FIG. 2. (Color online) Time evolution of the interaction energy Eint (Eq. (5)), the sum of the total kinetic energy and quantum 
energies (Eqs. (6) and (7)) and the dimensionless nonlinearity level, R (Eq. (16)). All energies are normalized with the particle 
number A^, as shown in Eqs. (10). Left panels: g — 1, Middle panels: g = 2, Right panels: g = 4. Upper panels: initial 
evolution, t < 5, Middle panels: intermediate stage, t < 40. Lower panels: latest stage, t < 4500 for g = 1, t < 2000 for g — 2, 
and t < 600 for g = 4. 



maximal is R ~ 0.2, as expected. An important point 
is that R is essentially less than unity (e.g. R ^ 0.2 for 
g 4 is the highest nonlinearity level at times 1.5 < 
t < 5) and decreases with time. This means that in an 
analytical approach to the problem, the interaction can 
be treated through a perturbation approach with respect 
of R. 



B. Evolution of the compressible and 
incompressible kinetic energy spectra 

The energy distribution between scales (energy spec- 
tra) for g — 1, 2, 4 and at different moments of time dur- 
ing free decay are shown in Fig. 3 for the compressible 
subsystem and in Fig. 4 for the incompressible subsys- 
tem. We conclude that: 

• At earlier stages (upper panels in both Figs. 3 and 
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4) the compressible and incompressible kinetic energies, 
being initially located at small k < 3Afc, propagate to- 
ward larger fc, still (almost) monotonically decaying to- 
ward large k. 

• The intermediate moments of time (around t ^ 3) 
will be discussed in detail in the next Sec. HID. Here 
we only note that the incompressible spectra are close to 
the Kolmogorov-Obukhov law El.i^{k) cx k~^/^ shown in 
the upper panels of Fig. 4. This behavior is related to 
the energy cascade. 

• At later times the energy begins to accumulate at 
large k and the energy spectra asymptotically approach 
the quasi-stationary state (practically coinciding e.g. for 
t ~ Ih and t = 40 within the error bars, related to the 
lack of statistical averaging). 

• The compressible energy spectra £'JJi„(A:) at this 
stage are close to the thermodynamic equilibrium with 
the energy equipartition between degrees of freedom. 
In 2D systems this gives (for one-dimensional spectra) 
^kin(^) ^ ^' shown in the lower panels in Fig. 3. This 
state will be discussed in more detail below in Sec. HIE. 

• The achieved (quasi-stationary) distribution of 
the incompressible kinetic energy El^:^^ (lower panels in 
Fig. 4) does not correspond to the thermodynamic equi- 



librium, as we can naively expect. This is related to the 
fact that at these times the system does not achieve full 
equilibrium. Indeed, the lower panels in Fig. 1 show that 
the incompressible kinetic energy continues to converge 
to the compressible energy. We interpret this stage as 
a kind of flux (not thermodynamic) equilibrium, when 
the spectrum E^^j^(k) is determined by the energy flux 
from the incompressible to the compressible subsystem. 
We think that at these times the energy exchange within 
these subsystems (that should lead to thermodynamic 
equilibrium) plays a subdominant role. 

• After a long time evolution, all the kinetic energy 
finally comes to comprise the compressible component. 
This is the full thermodynamic state. We show the 
long-time averaged compressible kinetic energy spectra 
-^kin(^)- expected, this is proportional to k. 

C. Evolution of the particle number spectra 

As we can see from Eq. (3b), the total particle num- 
ber is quadratic in \['(r) and therefore can be presented 
in the fc-representation by Eq. (10c), shown in Fig. 5 
for different time intervals and all values oi g — 1 , 2 
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FIG. 4. (Color online) Log-log plots of the spectra of incompressible kinetic energy at earlier (upper panels) and later (lower 
panels) moments of time. Left panels: g — 1, Middle panels: g — 2, Right panels: g ~ 4: 



and 4. As we see, initially, N{k) is localized in small-A:. 
Through the turbulent state, the spectrum asymptoti- 
cally approaches the form N{k) oc in the wavenum- 
ber range k < 27r/^. To rationalize this behavior, we note 
that in the full thermodynamic equilibrium this depen- 
dence should be related to E^[^{k) cx k. Indeed, in this 
state, the nonlinearity level Eq. (16) is small, meaning 
that the interaction energy can be neglected in compar- 
ison with the kinetic and quantum energy. In addition, 
the quantum energy which cannot be considered negligi- 
ble corresponds to zero-point motion and this term does 
not give rise to particle currents[7]. Thus, only the ki- 
netic energy spectrum is related to the particle spectrum. 
Moreover, almost all the kinetic energy becomes the com- 
pressible component and there are no vortices, and as a 
result has no singularities and becomes of order unity. 
Thus the kinetic energy density can be estimated using 
a wavenumbcr k as follows: 

^kin ^ f^in ^ Ip\VO\' ^ lpk'\9\' ^ \pe, (17) 



which becomes in fc-space using (10a) and (10c), 
\\w{k)\'^ ~ ifc2|^'(fc)|2, as a resuh, 

i?ki„(fc) ^ \k^N{k). (18) 

In this way, the two relations Ey^^{k) cx k and N{k) cx 
k~^ hold simultaneously. This relation is similar to the 
relation between enstrophy and kinetic energy E in 2D 
classical fluids. 



D. Observation of 2D-direct energy cascade 

Returning to the discussion of Figs. 1 and 3, we note 
the following: 

• The incompressible energy spectra at these times 
can be interpreted as closed to the classical Kolmogorov- 
Obukhov distribution E^^^^lk) cx shown in the up- 
per panels of Fig. 4. 

• Bearing in mind that during the system evolution 
the kinetic energy clearly propagates from small k (where 
it was initially located in £'J.i„(fc)) towards large k we con- 
clude that at some region of the system parameters we 
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FIG. 5. (Color online) Log-log plots of the particle number spectra at earlier (upper panels) and later (lower panels) moments 
of time. Left panels: g = 1, Middle panels: g = 2, Right panels: <? = 4 



can observe a 2D-direct energy cascade, which has previ- 
ously been observed only in 3D hydrodynamic systems. 

• Numerical analysis of the incompressible kinetic en- 
ergy flux supports the above explanation. Using Eq. (15), 
we calculate the incompressible kinetic energy flux. Aver- 
aged data for a short time interval r = 0.20 with five sets 
of data are shown in Fig. 6. The flux e'(fc) takes positive 
values for 3Ak <k< 2n/^ at least for 2.50 < t < 4.95. 
This strongly supports the occurrence of a 2D direct en- 
ergy cascade. Moreover, owing to the second term in 
Eq. (15), especially effective for high wavenumbers, £'(fc) 
tends to vanish at fcmax- 

• The Kolmogorov energy spectrum £'ki,j(fc) oc 

was previously observed in decaying 2D turbulence in the 
GPE [10]. However, this observation was interpreted in 
terms of the known 2D turbulence inverse energy cas- 
cade, predicted by Kraichnan [20]. We think that this 
interpretation is mistaken. In particular, if this inter- 
pretation is true, this spectrum has to be followed by a 
spectrum for a direct enstrophy cascade according to 
Kraichnan's scenario. However, the authors of Ref. [10] 
"observe a dependence in this range" , instead of k^^ 
dependence. 

• When the Kolmogorov-Obukhov spectrum is 



formed, a "two-dimensional" Richardson cascade can be 
seen. We can identify the cores of the quantized vor- 
tices by finding the phase defects of the wave function. 
Then, we choose the shortest intervortex pair length Ip 
for all vortices. In Fig. 7, the averaged vortex pair num- 
ber A/'pairCp) is shown and is proportional to /~", where n 
depends on g and 1.30 < n < 2.12. This power law sug- 
gests a self-similar spatial structure. For two-dimensional 
vortices, one of the most effective lengths, corresponding 
to the length of three-dimensional vortex ring, is the in- 
tervortex length. We conclude that Fig. 7 implies the 
following: First, vortex pairs with distances comparable 
to 27r/(3Afc) are created. Second, the vortex pair dis- 
tance progressively decreases. Finally, when the vortex 
pair distance becomes comparable to ^, the pair is anni- 
hilated. 



E. Frequency power spectra and types of motions 

1. Compresstble type of motion 

Important information about types of motion can be 
extracted from the frequency power spectrum, which is 
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FIG. 6. (Color online) Averaged incompressible kinetic energy flux in a short 
2.50 < t < 4.95. Left panels: g = 1, Middle panels: g = 2, Right panels: (7 = 4 



time interval r — 0.20 in the time range 
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FIG. 7. (Color online) Averaged vortex pair number as the function of intervortex length Ip. Left panels: g = 1, Middle panels: 
g — 2, Right panels: g — 4: 



the Fourier transform of the different-time pair correla- 
tion function of the motion amplitude. For example, if 
an amplitude A(t) oscillates with a particular frequency 
woi i-e. A{t) = Ao exp{iujQt) , then 

J(t) = {A{t + T)A*{t)) = Al exp(iwoT) , 



giving 

J{uj) = J J{t) exp(— iwt)c?t oc (5(w — wq) . 

In other words, pure periodic motion with frequency wq 
has a frequency power spectrum with a very intensive 
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peak of zero width at uj — luq. It is easy to check that 
decaying oscillations A{t) = Aq cxp(— 7< + iujot) have the 
power spectrum 



Bogoliubov [7] with the result: 



J{U!) OC 



7 



(W - Wo)2 + 72 ' 



(19) 



i.e. a peak at frequency ujq and width 7, inversely pro- 
portional to the life time of the motion. 

We observed exactly this kind of frequency power spec- 
tra for the compressible velocity component at later 
times, as shown in Fig. 8 for k — 3\/2 and k = 6\/2- 
As the system evolves to the thermodynamic equilibrium 
state, the power spectrum forms sharp peaks. This ten- 
dency is also seen in the frequency power spectrum of the 
wave function (not shown). We repeated this analy- 
sis for different k and in Fig. 9 we plot the position of 
the maxima of the frequency power spectra of compress- 
ible motion, averaged over a long time interval in the 
thermodynamic equilibrium state. Bearing in mind that 
for these periods the nonlinearity level R is very small 
(about 0.04 for g = 4 and even smaller for g = 2 and 1), 
so one can consider these fluctuations as almost linear 
perturbations on the background of the rest state. The 
Eigenfrequencies of these oscillations have been found by 



(20) 



plotted in Fig. 9. The excellent agreement between the 
theoretical and numerical results indicates that the ob- 
served thermodynamical fluctuation of the compressible 
velocity component is indeed Bogoliubov's elementary 
excitations. The relatively small, but finite width of the 
observed peaks characterizes the finiteness of the life time 
of these fluctuations, caused by interaction of the fluctu- 
ations with different fc- vectors. 



2. Incompressible type of motion 

What kind of frequency power spectrum can we expect 
for the incompressible type of motion? To answer this 
we note that an incompressible fluid exhibits vortex mo- 
tions in which the life time and turnover time are of the 
same order. This may mean that the position of the peak 
maximum ujq and its width 7 must be of the same order 
of magnitude. Moreover, more detailed analysis shows 
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that due to symmetry, = 0. In fully developed tur- 
bulence, the Kolmogorov-41 dimensional estimates gives 
T~^(fc) ~ e^/^k"^/^ for the inverse life-time of vortices, 
where e is the energy flux. 

Under a trivial consideration we can expect 7 ~ 
T~^{k). However, this will be true only in the reference 
system, co-moving with the small vortices. In the labo- 
ratory reference system, where the mean velocity of the 
fluid is absent, we have to account for the sweeping of 
small vortices in the velocity field V of the largest vor- 
tices. This gives a Doppler shift of the peak position 
equal to k-V . Noting that the turbulent velocity field is 
random, we average the shifted frequency peaks over the 
statistics of large-scale motions, in which positive and 
negative velocities have the same probability. As a re- 
sult we d o not shift the peak, but instead widen it to 
7(fc) ^ ./W^k. 

Qualitatively, this is what is seen in Fig. 10, where the 
frequency power spectrum of the incompressible effec- 
tive velocity component with k = 3\/2 and k = 6\/2 are 
shown. As expected, they do not form sharp peaks. Also 
the peak width for k = 3\/2 (upper panels) are smaller 
than for k = 6\/2, (1 ower panels). We consider these 
facts as additional support for the fact that incompress- 
ible motion demonstrates turbulent vortex behavior. 

As wc already demonstrated, El^^{t) approaches af- 
ter t 5, taking a long time, and all the vortices finally 
entirely disappear. As a result, the compressible effec- 
tive velocity component dominates in the thermodynamic 
equilibrium state. 

IV. CONCLUSION 

In this paper we predict and numerically demon- 
strate a two-dimensional direct energy cascade with the 
Kolmogorov-Obukhov's —5/3 law in two-dimensional 
quantum turbulence with the two-dimensional Gross- 
Pitaevskii model. The cascade is mainly caused by the 
compressibility of Bose-Einstein condensates which for- 
bids enstrophy conservation. A two-dimensional direct 
energy cascade has not previously been observed in typ- 
ical two-dimensional classical turbulence in the decaying 



case, which features direct enstrophy cascade and self- 
organization of the system, because of the conservation 
laws of both energy and enstrophy. 

In addition, we demonstrate that the turbulent state 
shifts to the thermodynamic equilibrium state full of el- 
ementary excitations without quantized vortices. The 
dispersion relation calculated from the frequency-power 
spectrum of the compressible velocity components shows 
that the equilibrium state is filled with Bogoliubov exci- 
tations. 

It is important to note that two-dimensional quantum 

turbulence is not restricted to theory. It is possible to 
design experiments of two-dimensional quantum turbu- 
lence in pseudo two-dimensional atomic Bose-Einstein 
condensates, superfluid "^He, and superfluid ^He. So 
far, almost all studies of quantum turbulence have con- 
centrated on the three-dimensional case, except for a 
few works [10, 26]. However, at least in atomic Bose- 
Einstein condensates, the compressibility causes a sig- 
nificant change from two-dimensional classical turbu- 
lence, and in a certain parameter range two-dimensional 
quantum turbulence differs from two-dimensional clas- 
sical turbulence. Thus, we are certain that the two- 
dimensional direct energy cascade in two-dimensional 
quantum turbulence provides a new phase of study of 
quantum turbulence. 
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